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ABSTRACT 

The methods involved in satellite tracking data gather- 
ing and the procedures followed in satellite orbit deter- 
mination from the tracking data are necessarily of interest 
both to the tracking systems engineer and the orbit analyst. 
It is the purpose of this paper to document for such inter- 
ests a description of the Goddard Range and Range Rate 
system's measurements^ limitations of these measure- 
ments, and effects of partial processing of the measure- 
ment data prior to final orbit computation. The description 
of measurement of satellite radial velocity, i.e. range rate, 
is given special attention since much controversy arises 
over its definition, measurement, and usage in orbit 
determination. 
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A REVIEW OF GODDARD RANGE AND 
RANGE RATE SYSTEM MEASUREMENTS AND 
DATA PROCESSING TECHNIQUES 

T. J. Grenchik 
B. H. Putney 

Goddard Space Flight Center 


1.0 INTRODUCTION 


Many of the shortcomings of orbit determination result from the separation and semi- 
independent operation of each of many specialists responsible for a portion of the overall orbit de- 
termination process. The separations may be classified loosely as tracking system designer, 
tracking system operator, tracking system engineer (monitor), data preprocessor (data compressor), 
and orbit analyst. The bridging and interchange of information between these disciplines would be 
the ideal means of error reduction in orbit determination and prediction. 

In this report the data preprocessing and the basic tracking system operation are described 
in simple terms for both the Goddard Range and Range Rate (GRARR) System and the Applications 
Technology Range and Range Rate System, an outgrowth of GRARR. In both systems the primary 
measurements are radial range and rate of change of radial range. In this report all significant 
system error sources are described and/or referenced. Additional tutorial references are in- 
cluded for those interested. Time -tagging errors are given the greatest attention because it is 
felt that this class of errors is one of the simplest to remove in the tracking data preprocessing 
(data editing and/or compression). Some representative values of time-tag errors are used to 
demonstrate their influence in orbit determination. 

In the time-tag error class, the uncertainty in WWV transmission time from the Colorado or 
Hawaii transmission sites to any STADAN tracking site is the dominant source of error for the 
range rate measurement. The uncertainty in the time of the reception at the station is approximately 
1 millisecond. If a constant radial acceleration of 400 m/sec 2 (worst case assumption) coincided 
with this uncertainty in time tag,, a range rate error of approximately 40 cm/ sec would result. It 
is estimated that equipment delays account for range rate time-tag errors of only 200 microseconds. 

In the range measurement, if a constant spacecraft radial velocity of 10 4 m/sec (worst case) 
is assumed to coincide with the 1 millisecond WWV uncertainty, a range error of 10 meters would 
result. For the range time tagging, there may be equipment time delay errors which are greater 
than 1 millisecond. In one specific case an error of 2 milliseconds was attributed to equipment 
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delays, Here, additional field measurement and laboratory simulation is called for to verify or 
disprove this magnitude of equipment delay. 

The range rate measurement has been described in more detail than the range measurement 
to emphasize its nature as an average measurement of range rate over a relatively long time 
interval. The range measurement can be categorized as an instantaneous measurement and does 
not require an averaging correction such as applied to the range rate measurement. 

Error sources such as speed of light uncertainties, ionospheric and tropospheric effects have 
not been considered in this report because they are beyond the specialities of the authors. 


1.1 PREVIOUS ANALYSES 

Many documents concerning Goddard Range and Range Rate (GRARR) System operation have 
been published. This paper repeats only that information presented previously which is necessary 
for the reader's overall grasp of system operation beginning with the recorded tracking data and 
ending just prior to orbit computation. In those sections dealing with processing of the GRARR 
measurement data for orbit computation little of what is included herein has been published pre- 
viously, This documentation of the data handling then, will aid the tracking system engineer in 
comprehending the importance of the various classes of errors in the orbit computation. 


1.2 TYPES OF GRARR MEASUREMENTS 

The Goddard Range and Range Rate System, as the name implies, measures range and the rate 
of change of range between the ground tracking station and spacecraft. These measurements are 
taken periodically (generally once per second) and labelled with the time of occurrence of each 
measurement. Hence they are time samples of the range and range rate. Along with the primary 
measurements of range and range rate, samples of antenna pointing angle are also taken although 
at present they are not generally used in orbit computation. A detailed description of the Goddard 
Range and Range Rate System, its conception, evolution, and performance is available in Refer- 
ence 1, but sufficient explanation of the measurements will be given in appropriate sections of this 
report for reader comprehension without delving into system detail. 

2.0 THE RANGE RATE MEASUREMENT * 

The range rate measurement is a measure of the apparent frequency shift of the ground and 
spacecraft transmitter frequency caused by the relative motion between the spacecraft and the 
tracking station. One may argue whether the Doppler effect is better described as an apparent 
shift in frequency or apparent shift in phase angle of the electromagnetic radiation. No difficulty 
arises from either consideration since phase and frequency are linearly related by: 

at = " (i). 
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where 0 is phase angle in radians and is radian frequency. Dr, Kruger formerly of Goddard Space 
Flight Center (GSFC) in Reference 6 amply treats the subject of Doppler frequency shift including 
relativistic effects. Here the subjec 4 - is treated more simply but the same results are attained, 

2.1 RANGE RATE AVERAGING ERRORS 

The most important point, which cannot be stressed too strongly, is that the measurement of 
the frequency suift requires a finite amount of time. In fact longer times are preferable since 
generally the frequency measure becomes better with the length of the phase count measurement. 

Of course, what must be considered with the longer measurement is any change in range rate dur- 
ing the measurement interval masked by the averaging effect. 

In the Goddard Range and Range Rate (GRARR) system, the average Doppler frequency, I (1 , 
is established by measuring the time interval required to count a predetermined number of cycles 
of the instantaneous Doppler frequency plus a fixed bias frequency. A bias frequency is summed 
with the Doppler frequency in order to sense the direction of the range rate in a simple manner. 

In the GRARR system & RR is measured where 


and S RR is the time interval in seconds required for counting a fixed number (n) of cycles of bias 
frequency plus Doppler frequency (f b +7 d ) . 8 RR is measured by means of a time interval counter 
which counts the number of cycles of a reference frequency occurring during the time interval 8 RR 
corresponding to N cycles of f b +T d . The readout of the time interval counter, c 0 , at the end of 
the measurement is related to S RR by: 


G o “ f r ^RR (3) 

where C 0 is the count of the time interval counter in cycles and f r is the reference frequency 
counted in cycles per second. N and the sample rate* are chosen such that 8 RR is between one half 
and one second. 

To estimate the change in instantaneous radial velocity that can occur during a one second 
interval, it is required to look at the second and higher order derivatives of radial range (d 2 R/dt 2 , 
d 3 R/dt 3 , etc.) under some assumed worse case conditions. In Reference 2, Section V, an elliptical 
orbit with an apogee of 60,000 nautical miles and a perigee of 150 nautical miles is assumed. At 
perigee the rate of change of radial velocity d 2 R/dt 2 can be a maximum of 435 m/sec 2 with ref- 
erence to a ground station almost underneath the point of perigee. Over a one second interval, if 
d 2 R/dt 2 were constant at 435 m/sec 2 the radial velocity dR/dt would change 435 m/sec. 

?Sample rate can be set at other than a one per second rate. See Reference 3. 
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For example, the counting of Doppler frequency (plus bias frequency) over a one second meas- 
urement interval centered at perigee would give an avei’age Doppler frequency of approximately 
zero cycles. The sign of the Doppler would be one polarity during the first half-second approaching 
perigee and change to the opposite polarity while receding from perigee. The average Doppler fre- 
quency, while correct, masks the instantaneous change in Doppler frequency. If the Doppler meas- 
urement were not made with the time interval centered at perigee, it is obvious how large differ- 
ences between the instantaneous Doppler and the averaged Doppler (represented by the count, C 0 ) 
would occur. Dr. Kruger in Reference 7 explores averaging errors for low circular orbits and 
suggests corrections to be made for all types of orbits. This correction requires estimates of 
higher order time derivatives of range available from nominal orbit predictions. This correction 
technique will be discussed in greater detail under the data handling section of this report. 

It should be noted that the value of d 2 R/dt 2 is a critical function of orbit geometry. For ex- 
ample, at the 60,000 nautical mile apogee of the previously assumed orbit, d 2 R/dt 2 is -1.85(10" 2 ) 
m/sec 2 (minimum value) so that a range rate change over a one second time interval in this case 
would amount to approximately 2 cm/sec. For the case of constant spacecraft radial velocity with 
respect to a ground station (approximated during certain trajectories) the length of the time inter- 
val during the Doppler plus bias measurement would have no effect on the computation of instan- 
taneous Doppler frequency. 


2.2 BASIC ELECTRONIC SCHEME AND THE DOPPLER EQUATION 

In Equation 2, Y d was defined as the average Doppler frequency without demonstrating its re- 
lationship to spacecraft average range rate,* In this section the relationship will be derived. Fig- 
ure 1 shows the Doppler extraction process as it occurs in the Motorola built GRARR system. A 
more detailed description is presented in Reference 2. This basic example is suitable for explana- 
tion of all GRARR system Doppler extraction processes. 

In the lower left hand corner of Figure 1 an oscillator of frequency f t /M is multiplied by m to 
arrive at a ground transmitter frequency of f t . At the spacecraft transponder an apparent fre- 
quency of Kf t is received because of the spacecraft motion relative to the ground transmitter. 
From classical physics the multiplier k for a fixed source and a moving observer (spacecraft) is 
given by: 

dR 


where K is a time varying dimensionless scalar, c is the speed of light in m/sec, and dR/dt is the 
instantaneous range rate of the spacecraft in m/sec. In Figure 1, v and dR/dt are used interchange- 
ably. The apparent frequency will be greater than the transmitted frequency, f t , under the conditions 

*In preceding paragraphs the measurements are termed radial range and radial velocity. Common usage also ascribes the simpler names 
range and range rate. 
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Figure 1— Doppler frequency extraction process, 

of decreasing separation (decreasing range and hence decreasing range rate). This apparent change 
in frequency is analogous to the higher pitched tone heard from the locomotive diesel horn upon ap- 
proach and the lower pitched tone upon recession. This physical property implies that dR/dt is 
negative for approach and K - (c - dR/dt )/c > 1 under these conditions, 

Upon reception at the spacecraft, the frequency Kf t is mixed to a lower frequency Mf 0 -Kf t and 
subsequently is used to vary the phase of the spacecraft transmitted frequency Nf 0 . Nf 0 is derived 
from the spacecraft oscillator f fl and the spacecraft frequency multiplier N. Thus, in the space- 
craft transmitted output, a downlink frequency Nf 0 is modulated by a frequency, Mf 0 ~Kf t , which 
contains the Doppler information on the uplink frequency f t . 

Fourier analysis of the spacecraft downlink transmission enables the separation of the modu- 
lated signal into three significant frequency components: 


Nf 0 (carrier frequency) 
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Nf 0 Mf 0 « Kf t 


(lower sideband) 


Nf 0 + Mf 0 - Kf t (upper sideband) 


Those are shown at the transponder antenna in Figure 1. Each of these components upon reception 
at the ground appears to be multiplied by K' because of the relative motion of spacecraft and 
ground station, K' is defined as 


K' 


c 

~W 
c + ^ 


( 5 ) 


* 


whore K' is a time varying dimensionless scalar, c is the speed of light in m/sec and dR/dt is the 

instantaneous range rate of the spacecraft in m/sec as defined previously. The derivation of K' 

* 

follows from the frequency cliange apparent to a fixed observer (ground station) from a moving 
transmitter (spacecraft). Recall that an apparent increase of frequency comes about from a dimin- 
ishing separation. dR/dt must be negative ruder these conditions and K' n e/(c + v) > 1, 

In the ground receiver in Figure 1, the carrier term K' Nf 0 along with a 60 MHz term which 
sets ihe 1 st if frequency is reconstructed by the Voltage Controlled Oscillator and the times N 
multiplier. The upper and lower sideband frequencies K' [Nf 0 ± (Mf 0 - K» t )] are filtered off by the 
loop filter. Before filtering the upper and lower sidebands (now at the 60 MHz if) are sent to a 
60 MHz mixer which removes the 60 MHz term, leaving only a frequency term k' Mf 0 - KK' f t . This 
frequency term along with a 10 MHz term are reconstructed by the 2 nd voltage controlled oscillator. 


The ground transmitter oscillator f t /M is mixed with the reconstru cted downlink carrier term 
K' f Q - 60/N MHz, the resulting frequency is multiplied by M, a bias frequency f b is applied in the 
next mixer and finally in the last mixer, all frequency terms are removed except f b + (kk' ”1) f t . 
The spacecraft oscillator terms (f Q ) have, been removed and the resulting expression (kk' - 1) f t 
when converted in terms of dR/dt equals 


f 


d 


dR , v 

dR ( r\ 

dt* * 

dF V * t7 

dR 

/ R\ 

c + dF 

* + c 


*1 


( 6 ) 


The time difference of r/c occurs because K affects the uplink at one instant and K' affects the down- 
link at time R/c later. This indicated quantity, f d , is a time varying function which is averaged by 
the measurement system over a time S RR . The average value of f d over the time $ RR becomes 
approximately, 


f 


d 


„ AH 
2 At 
AR 
c + At 


( 7 ) 


In review K, K' , and dR/dt were defined in terms of a decreasing range causing an increase in 
frequency, with dR/dt defined as negative under these conditions. By definition R decreases as the 
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spacecraft approaches the ground station. dR/dt must agree in sign convention with R, and is nega- 
tive for spacecraft approach. The quantity; 


(KK* - l) 


c 


2 


+ 


ar 

At 

AR 

At 


( 8 ) 


is greater than one and the Doppler frequency is greater than zero, the physical result expected 
from the Doppler effect. Combining Equations 7, 8, and 2, one obtains: 


JL 

'RR " „ AR' 

, ' 2 ST , 

f b + Sr f . 


C + 


At 


( 9 ) 


In terms of the, readout of the GRARR system range rate time interval counter (c 0 in Equation 3), 
Equation 9 becomes; 


C, 


Nf r 

„ AR 
“ 2 

f b + — AR f e 


c + 


At 


( 10 ) 


2.3 TABULATION OF DOPPLER MEASUREMENT ERRORS 

References 2, 3, and 4 consider the effects of various errors in the range rate measurement 
and place magnitudes on these measurement errors. Here the intent is not to set exact quantitative 
limits on these errors but to catalog the errors according to source. For range rate these three 
classes are: 

1. All sources that cause phase (frequency') errors in the f b + f d measurement. 

2. All sources that cause time interval measurement errors in C 0 . 

3. All sources that cause error in identifying the correct time of the sample measurement. 

3.0 DISCUSSION OF DOPPLER MEASUREMENT ERRORS 

3.1 PHASE ERROR IN DOPPLER PLUS BIAS FREQUENCY 
Oscillator bistabilities 

In Figure 1, there are shown a large number of oscillators: two voltage controlled oscillators 
in phase locked loops, a transmitter oscillator and multipliers, a spacecraft oscillator and 
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multipliers, receiver oscillators at 50, 60, and 10 MHz, and a receiver oscillator which places the 
bias frequency into the processing system. Each of these oscillators contribute phase instability 
to the extracted frequency f b + f d since in each conversion a minute phase instability is transferred 
from the imperfect '"'“Qillator to the Doppler frequency. Magnitudes of errors to be expected from 
each of these oscillators are listed in References 2, 3, and 4, It suffices here to state that the 
greatest source of error comes from the voltage controlled oscillators in the phase locked loops. 
This is a direct result of the requirement for a pulling range of several hundred kilohertz. Phase 
stability and the frequency change capability of a voltage controlled oscillator are in direct opposi- 
tion to one another; as one gets better, the other becomes worse. 

Axis Crossing Detectors 

Axis crossing detectors are used throughout the GRARR system to identify times of zero 
crossings (phase reference points). If a perfect axis crossing detector could be built, it would 
sense accurately the occurrence of any input signal’s magnitude becoming and passing through 
zero. (Limitations of real axis crossing detectors will be discussed under Time Interval Meas- 
urement Errors.) In system operation the input signal’s magnitude has been corrupted by noise so 
that the true occurrences of zero magnitude have been shifted randomly in time. In some instances 
under sufficiently strong noise conditions, an occurrence of zero magnitude will be deleted, that is, 
if one were counting zero crossings of a sine wave to measure frequency, the count would be low 
by the number of zero crossings suppressed by the noise (actually 1/2 of the zero crossings since 
frequency equals 1/2 of the number of positive and negative going zero crossings). 

The phase (frequency) errors caused by system input noise are irreducible unless improve- 
ments can be made in the overall receiving system’s noise content, or the spacecraft radiated 
power can be increased sufficiently to force the system noise into a less competitive position. If 
this were done similar improvements would be required in the system oscillator instabilities to 
prevent those sources of error from becoming dominant. 

At present all sources of error are approximately balanced. If new mission requirements 
should impose a large improvement factor in range rate error, all sources of error would require 
proportionate improvements. 


3.2 TIME INTERVAL MEASUREMENT ERRORS 

In the preceding sections it was shown that there are a number of sources of error which affect 
the frequency (phase of the Doppler plus bias frequency (f b + f d ). Iu this section those errors af- 
fecting the measure of the corrupted f b + f d frequency will be discussed. 

Once the f b + f d frequency has been detected in the axis crossing detector, the "noise-free" 
resultant is counted in a counter which has a preset count limit, N. At the beginning of the count 
of N a start reference is generated. At the end of the count of N a stop reference is generated. Be- 
tween the start and the stop reference a time interval is measured. For emphasis a preceding 
description is repeated here. This time interval is measured by counting cycles of a fixed reference 
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frequency (f r ) from the beginning of the interval until the end of the interval, & RR . The number of 
cycles counted multiplied by the time required for one cycle to occur provides a measure of the 
time interval needed to count N cycles of f b + f d , that is if 


RR f +7 
r b T r d 

the number of cycles counted (c 0 ) of the reference frequency would equal 


( 12 ) 


(13) 


C Q is the system measurement and it is this measurement which is used in orbit computation. Ex- 
pressing d RR in terms of c 0 and f r , and f d in terms of dR/dt and c, we have: 


Nf.. 




- 2 


f b + 


c + 


AR 

At 

AR 

At 


Start and Stop Pulse Jitter 


(14) 


A perfect axis crossing detector as des- 
cribed previously recognizes the occurrence 
of the passage of the input signal through zero 
magnitude. In reality an axis crossing de- 
tector works on the principle of a large ampli- 
fication of the input signal prior to a low level 
threshold detection of the amplified signal. 
Figure 2 illustrates this method of detection. 
Shifts in the threshold level and especially 
circuit unbalance in the positive and negative 
thresholds give the effect of shifting the posi- 
tive and negative slopes of the reconstructed 
"noise-free" output. Since the start and stop 
pulses originate from the 1 st and N th cycle of 
the reconstructed f b + f d frequency, the jitter 
in the reconstruction shows as jitter in the 
start and stop pulses required for the interval 
measurements. 

Quantizing Error 

The time interval measurement measures 
complete cycles of the reference frequency, 




NASA-GSFC-T&DS 

MISSION a TRAJECTORY analysis division 
BRANCH 551 DATE MAR, 5, |969_ 

nv T. GRENCHIK pi OT Nin 1153 


Figure 2— Axis crossing detector operation. 
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that is if the reference frequency f r is 10 MHz, tire maximum uncertainty in the time interval 
measurement is related to the time f r takes to complete a full cycle. The start pulse and the stop 
pulse are not coherent to the time base frequency f r so that a maximum error of O.lu second can 
occur at the start pulse and O.lu second at the stop pulse. 

These quantization errors can be reduced by increasing the frequency of the reference fre- 
quency. At present however, the quantization errors are adequately minute. (See References 2 , 

3, and 4.) 

f r Time Base Reference Short Term Instability 

f r , the 10 MHz reference frequency has an uncertainty associated with its phase (frequency) 
caused by the minute amounts of noise and spurious signal contained in its generation and ampli- 
fication. This is evidenced by lengthening or shortening of the cycles of 10 MHz signal which is 
counted for C Q . Limits on the uncertainty in f r are set by the mission requirements and the un- 
certainties generated by other error sources. 


3.3 TIME TAGGING ERRORS 


In Section 2.1, the averaging affect of the Doppler frequency measurement was noted. In this 
section similar types of errors will be described, all of which are present if the acceleration 
and/or higher order derivatives of range are not zero. Figure 3 shows graphically the significant 



Figure 3— Range rate measurement time tag. 
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times and time intervals associated with the range rate measurement. Definitions are as 
follows: 


t r - Time associated with range Rj 



Range at t r 

Spacecraft transponder delay at T R which is a function of uplink Doppler 
frequency 


T d - Data time as reported in GRARR format for measurement of range rate 

AT ft - Time interval following t d , after which range rate measurement actually 
starts. The reason for this delay will be explained later 

- Time of receipt at ground antenna of first cycle of spacecraft transmission 
used in Doppler frequency measurement 

AE^ - Time delay from antenna to Doppler extraction process. 1 st cycle of space- 
craft transmission reaches Doppler extractor at T + AE^ 

T rr - Middle of range rate measurement interval (spacecraft time) with which 
range rate measurement is associated in orbit computation 


T R - Time associated with range r 2 
R 2 - Range at l R 

- Spacecraft transponder delay at t which is a function of uplink Doppler 
frequency 

T 2 - Time of receipt at ground antenna of last cycle of spacecraft transmission 
used in Doppler frequency measurement 

S RR - Time interval measurement of n count of f b + f d 


Ar dR/dt (t r ) 


If T rr is not correctly determined the range rate measurement will have error since the acceleration 
and higher derivatives of the range are normally present. Some idea of magnitude of these time 
tagging errors, error in identifying the time t rr absolutely, are of interest. 

The times T x and T 2 in Figure 3 show the beginning and end points of the range rate counting 
interval S RR when referred to reception at the ground antenna. There is a finite time delay, AE^, in 
the equipment that extracts f b + f d , so that the counting process does not begin until Tj + AE^ , and 
the counting process does not end until T 2 + AE^ . The time on the ground must be referred to antenna 
input, actually to the Y-axis of the receiving antenna. (See Reference 8 for the timing reference 
point of the GRARR receiving system.) Hence, the times of the counting interval referred to the 
antennas input are T and T 2 . 
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t d is the reported data time of the measurement but there is a finite time delay AT^ after T D 
before the & RR measurement begins. This delay is caused by the first cycle of f b + f d clearing the 
N counter before proceeding with the measurement and counting to N + 1 cycles. 

In Figure 3, it can be seen that the cycle counting process is actually counting cycles of Dop- 
pler frequency which have occurred at the spacecraft prior to reception on tire ground. To relate 
the reception of the frequency at time to the time of the frequency at the spacecraft, T R , we 
have: 


t, • 


R. 


(15) 


Similarily we find: 


At 

t r 2 + T 


[ar (t r J 


R„ 


(16) 


The solutions for T RR follows from: 


T + T 
R 1 R 2 


RR 


T + T 1 

1 2 1 


RR 


2 4 1^ r i dt 


dR 


( T *,)] + AT [^( T *j]}‘ ^ K + M 


(17) 

(18) 


T, and T„ are related to T n and by: 

1 2 D 


T - T + AE- - AT- 

D 1 R R 


T, = T , + S RR 


(19) 

( 20 ) 


Equation 18 becomes upon elimination of T x and T 2 



At this point it is also necessary to insert a correction in Equation 21 for synchronization of 
the ground station time clock to the high frequency (HF) time transmission of WWV and WWVH. 
(See Reference 11.) For ease of synchronization, the ground station clocks are matched in time to 
the reception of the transmission and not to the time of the origination of the transmission.* Thus, 

*In the near future ground station clocks will be set on time, that is the clocks will include the time delay correction. When this occurs 
no WWV corrections will be necessary in the GRARlil data preprocessing. 
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for each station, its time clock is "slow" from tens to hundreds of milliseconds, depending on the 
time delay before reception of the time transmission. The data time t d used at the ground station 
is accordingly "slow" and Equation 21 corrected for this slowness is: 


RR 


+ WWV + AT,-, + 


‘D 


R 


'RR 

2 


AE/ 


1 (\ fdR 

~ 4 | Ar [HT 


( t O] +At [^( t O]}'^ (V* a ) 


( 22 ) 


where WWV is the time delay in seconds from the transmission of the WWV signal to the reception 
of that signal at the ground station. Each ground station has a measured value for WWV which is 
an average value of the time delay calculated from many observations. These time delay values 
are judged to be correct only within one millisecond because of the uncertainties in the height of 
the ionosphere, tilting of the ionosphere, etc. which cause variation in the mean path length traveled 
by the transmission. 

Recently, an experiment with a portable clock carried aboard the GEOS-B satellite has demon- 
strated a capability of synchronizing ground stations to within 0.1 millisecond, a factor of ten im- 
provement over the WWV transmission time uncertainty. This improvement was possible because 
of the precise definition of the GEOS-B orbit and the relatively smaller uncertainties in the time 
delay of the 136 MHz transmission carrier of the timing signals from the spacecraft to the ground 
station. 

These preliminary descriptions now allow the assignment of magnitude of time tag errors for 
T rr in Equation 22. WWV can be in error from 0.1 to 1 millisecond. The average value of AT^ is 
50 m seconds for the VHF measurement and 3 m seconds for the S-band measurement. AE^ is esti- 
mated to be 0.2 milliseconds (and will be discussed in detail in following paragraphs). At trans- 
ponder delay is never larger than 18 m seconds* so that 



is never larger than 9 m seconds. 

The term (r x + r 2 )/(2c) depends on the range to the spacecraft. For tracking a satellite orbit- 
ing the moon this term would be approximately 1.25 seconds. 

A review of the time tag errors shows that the order of importance of the errors is as follows: 


Maximum Error 
in Time Tag 

• lc ( R i + R 2 ) > ^ not considered in data 1,250, 000m sec 

processing 


Range Rage Error if 
d 2 R/dt 2 = 435 m/sec 2 

543 m/sec 


*VHF transponder delay is approximately 17 m seconds, S-band transponder delay is 

n 1 


approximately 3M seconds. 
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• WWV, uncertainty 

• AE^, if not considered in data 

pi'ocessing 

• AT^, if not considered in data 

processing 


Maximum Error 
in Time Tag 

100-1000m sec 


Range Rate Error if 
d 2 R/dt 2 = 435 m/se c 2 

.0435-. 435 m/sec 


200m sec 


.0870 m/sec 


3-50m sec 


.0013“, 021 m/sec 




i{Ar[dR/dt(T Ri )] 
•KArpR/dt (T Rj )]} 


if not considered in 
data processing 


9m sec 


,004 m/sec 


Consider a time tag error of 200m seconds such as AE^ . Again it is necessary to refer to the 
assumed worst case orbit condition with an acceleration, d 2 R./dt 2 , maximum of 435 m/sec 2 . Over 
a period of 200m seconds a constant acceleration of 435 meters per second squared would produce 
a range rate change of 0.087 meters per second, that is if a time tag error of 200m seconds coin- 
cided with a constant spacecraft acceleration of 435 m/sec 2 , the velocity measurement would be 
in error by 0.087 m/sec. 

The time tag error of AE^ is of special interest since next to the uncertainty in the WWV trans- 
mission time, it is the most significant possible source of error. Recall that AE^ is the time delay 
through the receiving system from the antenna input to the "n m counter input. In reexamination of 
Figure 1, nothing is shown which could cause significant delay through the system. Admittedly 
there is a physical separation of several hundred feet between the receive antenna and the "N" 
counter but this would only account for a time delay on the order of 1/3 of a microsecond. Yet 
time delays of approximately 200m seconds have been seen in one specific, case described in Ref- 
erence 9. In Figure 1 the phase locked loops themselves cause little or insignificant delay to the 
received signal, principally because the selectivity of the loop filter is small and phase change 
through the loop filter is correspondingly small. (See Reference 10 for description of the phase 
slope characteristics of phase locked loops.) What is not shown in Figure 1 are the 10 kHz band- 
pass filters placed ahead of the phase detectors in the phase locked loops so as to limit the noise 
bandwidth and present a reasonable signal to noise ratio input to the phase detectors. Gardner, 
et al., in Reference 10, describes the limitations of practical phase detectors which force additional 
filtering of the phase detector input. It is this additional filtering which introduces considerable 
delay through the phase locked loop. In any selective network, the phase angle of the input signal 
varies as a function of the input frequency, and from greater selectivity greater phase variation 
(delay) results. 

Corrections for this delay, if necessary, can only be obtained from measurements made on 
the phase characteristics of the selective circuit preceding the phase detector input. In the situa- 
tion described in Reference 9 where a consistent 200 microsecond delay existed in the range rate 
measurement time, it is certainly possible to presuppose a delay of this magnitude. No measure- 
ments to date have been made of any parts of the system delay affecting the range rate measurement 
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but future investigation will be made to assure that this source of error lias not become 
predominant. 

Actual corrections made in the data processing for time tag errors are described under the 
data processing section. This section is only meant to order the importance of the errors and con 
sider their effect on the range rate measurement error. If range rate measurement errors for 
frequency (f b + f d ) and time interval S RR are greatly reduced, then time tag errors will require 
similar reductions. At present, the uncertainty in the WWV transmission time is the greatest un- 
known source of error in time tagging. 


4.° THE RANGE MEASUREMENT 

The range measurement in a measure of the time interval between the transmission of a time 
marker or reference and the reception of that time marker returned from the cooperative target. 

The time marker in the case of the GRARR system is a mark or reference on the phase of the 
ground oscillator, f t , Upon reception of that phase marking the intervening time interval is meas- 
ured and represents the time delay for the reference mark to travel to and from the cooperative 
target. Given a constant velocity, c, for the travel of the reference mark to and from the target, then: 

2 R s c § R (23) 

where R is range in meters, c is the constant speed of light in meters per second, and S R is the 
time interval in seconds.* 

Within the cooperative target the phase mark on f t in the uplink is transferred to a phase mark 
on the spacecraft frequency Nf 0 in the downlink. Associated with this action performed on the space- 
craft is a delay which prevents the turnaround of the phase marking from occurring in zero time. 


4.1 TRANSPONDER DELAY 

The required selectivity built into the spacecraft transponder imposes a delay on the passage 
of the time reference mark through the transponder. Variations of input frequency to the trans- 
ponder caused by the Doppler shift impose time delay variation of the time reference mark through 
the transponder. For a discussion of these points, see Reference 12. Hence Equation 23 must be 
modified to account for the time delay through the transponder, a function of the Doppler shifted 
uplink frequency, Kf t : 

2E = 4„-Ar[f(Kf t )]} (24) 


•No attempt is made here to consider the bending of the path of travel in the earth’s atmosphere, which makes the line of sight distance 
shorter than R, the measured distance. 
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Figure 4— Range measurement-. 


where At is the time delay of the time refer- 
ence mark through the transponder, For ex- 
actness however, Equation 24, should contain 
additional terms. These terms have been neg- 
lected as insignificant, but it is of interest to 
examine their magnitudes. 

Figure 4 shows the important events in 
the range measurement sequence. b Rf the time 
interval measurement of the system is: 


( I'r j “ T[ ) + At ♦ (T 2 -T Rj ) 


(25) 


Rj + R 2 


+ At 


(26) 


A simplification has been made by omitting the 
ambiguities in the measurement. This will be 
defined and covered in a later section. Just as 
for the range rate measurement, the change in 
range from the beginning to the end of the trans 
ponder delay should be considered to determine 
its effect upon the measurement. r x and r 2 can 
be expressed in terms of r m and t m by a Taylor 
series expansion of R in the neighborhood of 
the point t “ t m 


R(t) 


"(T„) * af (T u )[t-T„] + ■§ ( T J [ t -T M ]» + 


cP_R 

dt- 


(27) 


Use of Equation 27 for the purpose here is approximate for 2 reasons, one, that by truncating 
the series with a finite number of terms one neglects the contribution of higher order derivatives 
of r, and two, considering the neighborhood of the point t m to include the transponder delay Ar is 
not completely rigorous. In defense of the approximation in Equation 27, because of orbit dynamics 
it is very unlikely that a large number of higher order derivatives of r are significant over a short 
period of time. Secondly, the transponder delay, At , is never larger than 18/x seconds. 

To evaluate Equation 27 for R = R x and R = R 2 , insert the appropriate values for the time and 
range values into Equation 27 


r ( t m) + af ( T «) [T a - T„ ] + \ ( T „ ) [t, - t m] 


2 + 


(28a) 
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* t ft 


(28b) 


R 


2 


K(T«) * § (**)[*! T„] 


l ‘HR 
5 dt» 


(% l ) [^2 



* 


But 


ancl 


T, " T„ 



T 


2 


T 


Sr 

T 


(29a) 


(29b) 


so that 


R, * R s 


2R(T„ ) 


iHE fT \(££\* 

dt 2 \ l Ul\ 2 ) 


(30) 


and all terms except even order derivatives drop out because of the difference in sign on St/2, The 
assumed worst case value of d 2 R./dt 2 is 435 meters per second squared. Using this value for the 
second term in Equation 30, one obtains: 


d 2 R 
dt 2 

U 



,-V w.j.ue 


3,5 (l0“ 8 ) meters 


( 31 ) 


It is readily agreed that this value ir negligible when compared to the system accuracy re- 
quirements. Thus Equation 26 becomes: 


®(T m ) 


+ At 


(32) 


4.2 TABULATION OF RANGE MEASUREMENT ERRORS 

References 2, 3, and 4 consider the effects of various errors in the range measurement and 
place magnitudes on their measurement errors. Here as in the range rate description the intent 
is to catalog the errors according to source. For range the three classes of error are: 

• All sources that cause phase (frequency) errors in the phase marking of f t and Nf 0 . 

• All sources that cause time interval measurement errors for 8 R . 

• All sources that cause error in identifying the correct time of the sample measurement. 
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The first two sources of error are adequately covered in References 2, 3, and 4 and no further 
discussion will be given here, The 3rd source will be considered further. 


5,0 D ISCUSSION OF TIME TAG ERRORS IN RANGE MEASUREMENT 

Figure 5 illustrates the significant times and time intervals for the range measurement. Ref- 
erence 8 describes the method used in removing equipment delay times in the measurement of <> R , 
but Reference 8 docs not describe how this "zero set" of range affects the range data time tagging. 
Figure 5 attempts to show the displacement of the measurement from the time tag. Definitions are 
as follows; 

t d - Data time as reported for range measurement in GRARR format. 

AT r - Delay caused by filtering of phase marking tones after generation in time code generator. 

AT r - Time interval following T D + At f , after which range measurement actually starts. 

AE„ - Time delay from origination point of phase marking to transmission from antenna. 

At - Time interval for range to and from collimation tower (used in zero set) to system 
antennas. 

;> R - System measurement time for zero set procedures (precalculated). 

AE R ^ - Time delay from reception of spacecraft signal at ground antenna to completion of 
range measurement, o R . 

Reference 5 describes the mechanization of the phase marking procedure, In the GRARR 
Motorola system low frequency tones, 8 Hz, 32 Hz, 160 Hz, 800 Hz, 4000 Hz, 20,000 Hz, and 500,000 
Hz, are used to phase mark the transmit frequency. * The generation of these tones in the GRARR 
time code generator results in square wave signals which require filtering to remove the undesired 
higher harmonics. It is this filter delay which is described as AT F . The phase marking tones in 
square wave form are in synchronization with the data time T D since they are all generated in the 
same circuitry, the GRARR time code generator. After filtering the sinusoidal phase marking 
tones are delayed with respect to t d . 

AT R is a delay incurred in the start pulse generator which produces the signal to begin the 
measurement. The delay, AT R , is 25, 30, or 31/x seconds, depending upon whether the 500 kHz, 

100 kHz, or 20 kHz is the highest frequency used in modulating (reference marking) the uplink 
frequency, f t . 

AE Rj and AE Rj are all other delays occurring in the system between the equipment and the 
antennas. How these delays, and the others mentioned above, affect the time tagging will be explained 


♦A current modification program will also allow modulation with a psuedo random code for ambiguity resolution. 
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6_ MEASUREMENT WHICH MUST BE SET EQUAL TO 
‘ H AT FOR ZERO SET PROCEDURE 


AT_ AT, 


R 


AE n 


TIME INTERVAL FOR RANGE FROM COLLIMATION TOWER 
TO ANTENNAS 


AT 


H 


AE- 




DATA 

TIME 

T„ 


REFERENCE 

TIME 

OCCURS IN 
EQUIPMENT 


'REFERENCE 
TIME OCCURS 
AT 

TRANSMIT 

ANTENNA 


'REFERENCE 
RETURN AT 
EQUIPMENT 

REFERENCE 
RETURN AT 
RECEIVE ANTENNA 
FROM COLLIMATION TOWER 


8 MEASUREMENT SET EQUAL TO AT ABOVE 
TO ZERO SET SYSTEM 


ATr 


AT r AE +AE, 




AE = 


A Ed 


t: 


DATA 

TIME 

Tn 


ZERO SET 
REFERENCE 
TIME HAS 
BEEN DELAYED 
BY AE + AE, 


‘ZERO SET 
REFERENCE TIME 
AT ANTENNA 
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REFERENCE 
RETURN AT 
RECEIVE ANTENNA 
FROM COLLIMATION TOWER 


'REFERENCE 
RETURN AT 
EQUIPMENT 


Figure 5— Displacement of reference time from reference time with zero-set procedure. 

by reference to Figure 5. If the method of GRARR system zero-set is understood, the time tag 
errors are readily apparent. The measurement desired, S R , is the 2 way travel from transmit 
antenna to target and back to the receive antenna. (Recall that the Y axis of the receive antenna in 
the reference point of the system.) 

In the top figure of Figure 5, the data time t d occurs, and some time later after AT F and at r , 
the reference time is reached. This is the time at which the measurement S R begins. Some time 
later (after AE Rj ) the reference marked transmission leaves the transmit antenna, is turned around 
by a target on a collimation tower, and is returned to the receive antenna. This 2 way distance to 
the collimation tower (used in zero setting) is known very accurately by survey measurement. The 
returned transmission then passes through the equipment, (delay time AE R2 ) and completes the S R 
measurement. 
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Iii the lower figure of Figure 5, the <* R measurement which includes Ae r and AE r must be set 
equal to the accurately determined 2 way delay to the collimation tower in order to zero set the 
system. This is done as follows: The data time T d and the reference return at the equipment are 
fixed and unchangable. The only adjustment that can be made is adjustment of the reference time 
within the system. This is delayed a time, AE R + AE R ^ , so that the adjusted o R measureir^v now 
agrees exactly with the survey measurement to and from the collimation tower. 

It follows then that the time tag for the range measurement (spacecraft time) should be: 

T n " T„ * AT r + AT, * AE R| t AE Ri t ^ (33) 

Ambiguity terms have been omitted for clarity. Equation 33 requires adjustment of station data time 
T d to WWV time just as described under time tagging for the range rate measurement 


fc R 

T r = T D + WWV + AT p + AT R + AE Rj + AE R<> + -j (34) 

Except for the term At r , the proven magnitudes of At p , Ae r , and Ae r are unknown at this time. In 
Reference 9, a delay of 2.1 milliseconds was attributed to a time tag error in the GRARR system 
at Rosman, North Carolina. Verification of the time delays, AT f , Ae r , and Ae r remain for future 
work, 

A review of the time tag errors for the range measurement show the order of importance for 
the errors to be as follows: 



Maximum Error 
in Time Tag 

Maximum Range Error 
(dR/dt * 10 4 m/sec) 

S r / 2, if not considered in 
data processing 

1,250,000/x sec 
(one-way delay to moon) 

12,500 m 

AT*. + at r +ae Ri +ae Rz , if not 
considered in data 

2000^* sec 

20 m 

processing 



WWV uncertainty 

100-1000/u. sec 

1-10 m 


In the case of range, it appears that the WWV uncertainty is not the predominant error in the 
time tagging. Further effort will be expended in field measurement of the system delays. 


‘Value of time tag error taken to be equal to Roman bias (Reference 9). 
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6 *° PREPROCESSING of range, range rate, and angle tracking dat a 

6.1 INTRODUCTION 

Up to this point the discussion has been restricted to the field station measurement or the ac- 
quisition of the raw (uncorrected) range, range rate, and angle tracking data. In this portion of 
the document the viewpoint will be slanted toward the preprocessing of this raw data prior to its 
use in orbit computation. The process described herein includes the raw data reading, raw data 
smoothing, and conversion from range and range rate data in time intervals (seconds) to range and 
range rate values in meters and meters per second. The preprocessing methods do not include 
corrections for refraction. 


6.2 GENERAL DESCRIPTION AND PURPOSE OF PREPROCESSING 

The preprocessing program performs the following general functions. 

• It smooths by least squares fit to selectable order (up to 6th degree) Chebyshev polynomials 
large (up to 400 points) groups of raw data. The purpose of this large scale smoothing is to 
identify "bad" data points and to provide an estimate of data scatter (standard error of fit). 
Bad data points are identified by rejecting data points whose residuals exceed some factor 
of the standard error of fit, Iteration and rejection is continued until the percent change in 
the standard error reaches some specified number. 

• It computes processed data points for orbit determination using the following techniques: 
from the large groups of data smoothed, small groups of data points are selected. 

A single point is computed from the fitting polynomial at the center of the small group. 

These points are then converted to output (processed) points as described in the succeed- 
ing sections of this report. 


6„3 POLYNOMIAL FITTING TO RAW DATA 

The raw data measurements, in units of fx sec, are functions of time and are fit to polynomial 
expansions in time. A least squares solution using Chebyshev polynomials is performed on 
blocks of data points (nominally 96). 

6.3.1 Least Squares Solution 

Given is a set of n + 1 discrete data points |f(x 0 ), • • • , f(x n )}. The approximating function 
y(x) which approximates (fits) these data points is 


f ( x ) ~ s y( x ) 
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where s Q (x), • • • , <p n (x) are n + 1 appropriately chosen functions. Define the residual R(x) 


R(x) 


n 

f ( x > - ^ 
k=0 


f(x) - y(x) 


2 

= minimum , 

This is equivalent to saying 


for r = 0, • • • , n 

or 

for r = 0, • • • , n 

or 

for r = 0, • • • , n . The final equation describes n + 1 simultaneous equations in the n + 1 unknown 
parameters a 0 , a , • • • , a n . These equations are called the normal equations of the process. (See 
Ch. 7, Reference 14, Hildebrand.) 

Written in matrix notation we have de = F and d is a square symmetric matrix with the general 
element 

N 

d st * ( X l) 

4= 0 

s , t = o, » • • , n and F is a column vector with 

N 

f s = 2Z ^ 

4=0 

s = 0, 1, • • • ,n. Vector e is the solution, (e = D" 1 f). (See Reference 15, Berztiss.) 


<9 


-A 2 


\4=o L k “° 


= 0 


11 n 

^ h )- 7"a k 0 k (\fc) 
4=o 


k=0 


n N 


k=0 4=o 


and the least squares property is defined as 


11 n 

a 2Z f fa)~ 2T. a k^( x 4) 

4=o 4=o - 


k= 0 
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6.3.2 Choice of Chebyshev Polynomials as the Data Approximating Function 

Most often it is asked "Why not use the monomials (Taylor series expansion about x = 0)?" 
That is, let 0 O (x) = 1, 


(x) = x , ’ * * i <£ n (x) " x n 


and 


f(x) = a Q + a x x + • • * + a n x n , 

As described in the referenced Berztiss paper, the matrix D in the case of the monomial ex- 
pansion has a much greater tendency to be badly ill-conditioned than an expansion in terms of 
orthogonal polynomials. 

So one turns to orthogonal polynomials, which are defined as satisfying the condition 

f 0i ( x )0k ( x ) dx “ 0 , 


for the set of functions $ 0 (x), • • • , <£ n (x) , over the range a < x < b. The Chebyshev polynomials, T. (x), 
have the orthogonality property 


f 1 T r (x)T s (x) 

___ — dx = 0 (r / S') 

J - 1 y 1 - x 2 

which makes them orthogonal with respect to the weight factor l/fi - x 2 . 

In the special case when x 1 , • • • . x N are the n zeros of T n (x) , the off diagonal terms of the 
matrix D vanish and the Chebyshev polynomials are truly orthogonal. It ’ i a desirable property 
of orthogonal polynomials, that the d matrix becomes a diagonal matrix. Although the zeros of the 
T n (x) polynomial are never available, the resulting D matrix from the Chebyshev polynomial so- 
lution has comparatively small off diagonal terms and makes the solution well defined. According 
to Berztiss, Reference 15, the use of Chebyshev polynomials compared to monomials improved the 
accuracy for a 10th degree polynomial by approximately four significant figures using a single 
precision arithmetic in a computer. This difference between the Chebyshev and the monomial ac- 
curacy could have been made arbitrarily small if arithmetic of sufficiently high precision (i.e. 
double precision) had been used, but it does illustrate the advantages of minimizing the values of 
the off diagonal terms in the D matrix. 
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The Chebyshev polynomials used in the preprocessing are: 

n 

y (x) Mi (x) 

kc' 0 

where 

x = time normalized to [-1, 1] interval 
a k = the desired coefficients 
n = the degree of the polynomial 
T\ (x) = the Chebyshev polynomials 

where 

T 0 (x) - 1 

MX) s X 

and W*) a 2xT r (x) - MCx) 


6.4 ORB 1 TAPE 

Some of the difficulties inherent in the preprocessing of raw data stem from insufficient data 
available from raw data sent from the field stations. For example the reported range measurement 
may be ambiguous in that from the data message one cannot determine if the range is 18,750 kilom- 
eters or multiplies of 18,750 kilometers. This difficulty is circumvented by the prior preparation 
of an ORB 1 tape, which is used in conjunction with the data preprocessing. The ORB 1 tape is a 
record which contains an ephemeris of position and velocity of a satellite with reference to the 
center of the earth. 

There are several Orbit Generator Routines in the final orbit determination program that have 
as a by-product an ORB 1 tape. The Cowell Integration Orbit Generator Routine is the most fre- 
quently used at the present time. Needed as input to these routines is a set of orbital elements (co- 
ordinates that can be converted to position and velocity at that time) and a time to be associated 
with this set of elements. When a satellite is first launched, the elements used to create the first 
ORB 1 tape, are the predicted elements. When the first data comes in, it is in conjunction with 
these initial elements that the observations are differentially corrected in the final orbit determina- 
tion program and a new set of orbital elements are computed. With this new set of elements a new 
ORB 1 tape can be generated. For most satellites, after a satellite’s orbit has been fairly well de- 
termined, it is necessary to update the orbit about once a week. Regardless of the means of pre- 
paration, it is this ORB .1 tape which is used in' raw data preprocessing. Its specific usage in pre- 
processing will become more evident in later sections of this report. 
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6.5 RANGE MEASUREMENT DESCRIPTION 


The GRARR range measurement is a count of a reference frequency (100 MHz) from a start 
time (signal transmission) until a stop time (signal reception). Knowing the reference frequency 
this count can easily be converted into a time interval measurement. Using the speed of light one 
can convert the time interval measurement into a distance measurement. 

6.5.1 Raw Data Format 

The preprocessing program accepts the raw data in the formats governed by the Stadan En- 
gineering Division, Code 570. A description of the data formats is best available in a GSFC X Doc- 
ument to be published in April 1969, entitled, "DATA FORMATS OF THE GODDARD RANGE AND 
RANGE RATE SYSTEM AND THE APPLICATIONS TECHNOLOGY SATELLITE RANGE AND 
RANGE RATE SYSTEM," by DAVID ZILLIG, Code 571. 

6.5.2 Range Measurement Ambiguity* 

The time interval counter for the range measurement has a maximum value which is deter- 
mined by the lowest frequency sidetone (usually 8 Hz). The 100 megahertz counter is started on an 
even increment of time (usually once per sec- 
ond). The event that stops the counter is the 
return of a positive going zero crossing of the 
8 Hz modulating waveform upon return from 
the target. The time interval between two cy- 
cles of the 8 Hz sinusoid is .125 seconds. 

If a satellite’s range requires more than 
.125 seconds for the wave to be transmitted and 
received (a range of about 18,500 km) we say 
it is out of the first range gate and depending 
on which range gate it is in, a fixed quantity of 
time must be added to the measurement. 

To illustrate, let a time mark (8 Hz posi- 
tive going zero crossing) be sent at T D where 
T d is the time associated with the raw data. 

Figure 6 illustrates how ambiguity occurs. Let A indicate the positions of positive going zero 
crossings of the 8 Hz modulation on the transmitted carrier wave as the range measurement is 
occurring. The event that stops the range time interval counter is the reception of number "6" 
positive going zero crossing. Number "1" 8 Hz positive going zero crossing is beginning trans- 
mission at T d . The time measured S R is the time required for number 6 to complete its trip. 


REFERENCE POINT LEFT GROUND 
6 5 4 3 2 

J L 


,5 (COUNTER 
r R READING) 




TIME- 


0 6 5 4 3 2 1 

REFERENCE POINT RETURNED TO GROUND 


RANGE 



TRANSMISSION 


RECEPTION 


POSITION OF REFERENCE POINT SVJTH 
RESPECT TO SATELLITE AT Td TIME 

Figure 6— Ambiguity illustrations. 


*The ATSR system requires ambiguity resolution. 


In the near future all GRARR systems will have self-contained ambiguity resolution, 
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Therefore, the desired transmitted time is when number "6" left the ground and to co ^pute the 
proper range the number of gates or ambiguities must be added to the computed range. For this 
example it would be necessary to add .625 seconds (5 x .125) to the range time interval counter 
reading. 

6,5.3 Conversion of Range to Metric Quantity 

Let h R be the two way delay time (range measurement) in seconds, Ar be the transponder delay 
and c is the speed of light. From distance equals velocity multiplied by time one obtains: 

runge b k t ^-(b R -Ar) 


If an ambiguity is present, this equation must be modified to' 

R = § (& r -At +n a AA) 


where 

At b transponder delay 
N a - the ambiguity number 
AA = size of the ambiguity interval in time 
S R = two way delay time 

To determine n a use the ORB 1 tape described previously to compute range at t d where t d is the 
data time. Let 

R c ~ range computed from ORB 1 tape 

R i “ f S R , range computed from measurement ignoring ambiguity 
R a = the range associated with AA time -(gate size) 

Then let N a be the integer closest to X A when, 

x a r a = R c “ Rt or 
R c - R x 
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6.5.4 Range Measurement Time Tagging 

The data comes in tagged with a raw data time t d . To associate the proper data time with the 
measurement it is necessary to know the approximate orbit. The ORB 1 tape is used for this 
purpose. From the ORB 1 tape, which contains an ephemeris of geocentric position and velocity, 
a theoretical range with respect to a particular station can be computed. The desired time to be 
associated with the theoretical range is satellite time, that is the time of the range measurement 
at the satellite and not the time at the ground station. The data however, is time tagged at the 
ground, introducing errors because of the time required for the signal to reach the ground, the time 
required for the signal to go through the transponder, (transponder delay), the WWV time delay 
(defined earlier) as well as ground equipment delays. 


In order to determine the midpoint time of the entire measuring interval (from the time the 
signal leaves the ground to the time it is received) we use the ORB 1 tape to determine the approxi 
mate range. The transmitted signal time is unknown. The received signal time can be computed 
using the raw data tag, the measurement and the WWV correction. From the ORB 1 tape, the ap- 
proximate range at the received time can then be computed. Since r = ct and t = R./c , the mid- 
point time can be estimated, and the approximate range at the midpoint can be recomputed. This 
value is used to compute the midpoint time. The time is then corrected for the transponder delay. 


Let 


T d = raw data time 
R = the time measurement 
N a = the ambiguity number 
AA - the size of the ambiguity in time 
T r ta the proper data time tag 
WWV - WWV delay 
c a speed of light 

R c (t) =s range computed from the ORB 1 
tape at time t . 



s = T D + S R + wm T Jtli = t 2 - ±r c (T) 


As a first approximation use t 2 as an estimate for satellite time to compute r c at t r . 
As a second approximation 



This is sufficiently accurate at this time, and all that is necessary is to subtract for the transponder 
delay of the satellite. The transponder delay is a function of the range rate of the satellite. 
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Let A? (<1R( t ) dt) be the transponder delay as a function of range rate (JR <it computed from the 
ORB 1 tape at time t . 

Then 



is the time tag to be associated with the range measurement. 

6.6 RANGE RATE MEASUREMENT DESCRIPTION 

The range rate measurement is a count of a reference frequency (10 MHz)* starting with the 
first of n cycles of f b + f (1 (bias + doppler frequencies) and stopping upon the completion of the N th 
cycle. By knowing the reference frequency, one can interpret the measurement as a time interval 
measurement, a time interval required to count to a known value N of cycles of doppler and bias 
frequency. 

to 

6.6.1 Range Rate Formula Derivation and Averaging Effect Correction 

Let h m be the time it takes to measure N cycles of doppler frequency, f d , plus bias frequency, 
f b , in the range rate counter. 

N ' ( f d + f b) f, RB 

Since the measuring over n cycles takes a finite amount of time, the measurement is an average 
range difference. It is necessary to make a correction to the average range difference, using the 
ORB 1 tape, in order to make the measurement an instantaneous range difference. This is nec- 
essary since the orbit determination requires an instantaneous measurement, 

Let 

AR (b RR ) « the change in range during the ground time interval f? RR 
c = the speed of light 
f t s the uplink frequency 
f d = the 2 way doppler frequency 
f b * the bias frequency 

*For ATSR range rate, the reference frequency is 100 MHz. 
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which is the wave length of the transmitted signal multiplied by the average two way doppler fre- 
quency divided by 2, multiplied by the measurement time. The negative rlgn if necessary to main- 
tain the convention of negative doppler frequency corresponding to a positive change in range, 

Since 




HR 


( <l fc RR N " ' RR 


and "R(’ kb ) becomes 


■‘ r (‘rr) 2^( N - f u ; ’rr) 

The time to be associated with this change 
in range needs to be satellite time. 



m - 
ar 


1 a;, r 


liO. 




where Ac R is the change in two way delay time due to the change in range, Then for satellite time 


Ac, 


At 


RR 


At 


AR(At) 


RR 


s rr * 5T, (N'f b ‘ R *) 

since At at the satellite corresponds to & RR on the ground. Then: 

2t; ( n ~ f b ^rr ) 

^RR + 2f; (N" 'i^Rr) 


At 


which is the average range difference for that time interval. 
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To correct for the fact that this is an aver- 
age range difference the ORB 1 tape is used, The- 
oretical ranges are computed at t r (Rj) , and at 
t r 2 ( r 2 )> anci tllc theoretical range rate <1R <lt at 

(v^)/ 2 - 

The correction is then 


f 


corr 


dR / R 2 R i \ 
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and 


dR ARC At) 

H7 * vorr 


is the corrected instantaneous range rate. 


6,6,2 Range Rate Time Tagging 


The time associated with the range rate measurement is a reference time that needs to be ad- 
justed to satellite time, In essence the midpoint of the measuring interval is to be computed. Since 
the satellite is moving during the ground station measurement of the N cycles of f b + f d , It is nec- 
essary to determine where the satellite was at two points, the beginning and the end of the measur- 
ing interval, At these two points it is necessary to determine the satellite time, The procedure 
is the same as for the range measurement but two points in time are needed for the range rate 
measurement, 

Let 


RR 


RR 


t d a raw data time 

« the time measurement 
a the proper data time tag 
WWV = wwv delay 
c *8 speed of light 

R t . (t), dR(t)/dt a range, range rate computed from the ORB 1 tape at time t. 

T 2 = T n + + WWV 

Using the same procedure as uSed in the 
range computation and letting 

dR 



t 2 .. 2 - i (a£ (t,.,,)) 


¥ 
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Substituting T d for T 2 


> 


Then 


* 



T 


RR 




1 


Note that minor delays through the equipment are ignored. 


6.7 DATA OUTPUT FOR THE DIFFERENTIAL CORRECTION PROGRAM 

The actual procedure used to obtain representative points for the Differential Correction 
Program (final orbit determination) is as follows: 

c Smooth the raw data as a function of raw data time (equally spaced intervals) and obtain 
coefficients for the representative Chebyshev polynomials. Use all the data points and ob- 
tain a fit. Reject data points outside a predetermined number of sigma, (root mean square 
of the residuals) away from the curve. Compute new coefficients and continue in this process 
until sigma remains relatively constant or no new points are thrown out. 

• Choose output times and compute corresponding measurements from the polynomials. 

• Convert to range and range rate measurements. 

e Make time corrections. 

7.0 CONCLUSION 

It has been the purpose of this paper to shed some light on the range and range rate measure- 
ments and how they are processed prior to orbit computation. In particular, most errors of the 
system are tabulated and ranked. When correction is warranted, a method of correction is indi- 
cated. It is the hope of the authors that this description and partial analysis will prove useful to 
persons interested in the GRARR and ATSR tracking systems. 

Acknowledgment and thanks is made for the helpful suggestions in formula derivations and time 
corrections to Dr. Joseph Siry and Mr. David Stewart. 
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Appendix A 

Past Preprocessing Methods 

A.l TIME TAG PLACED ON RANGE AND RANGE RATE DATA 
ON THE IBM 7094 MACHINE 

As described earlier in this document, the time tag placed on the raw data at the tracking 
station is the same for range and range rate measurements, although the time tag should be in- 
terpreted differently. The time that this raw data time tag (t d ) represents can be seen from the 
diagrams. 




The transponder delay (At) in the IBM 7094 program was taken as a constant. The value 
chosen was the value at zero doppler frequency. The time tagged to both the range and range rate 
data by the preprocessor was: 


( § r“ A t ) 

T = T d -t- WWV + j 


expressed in compatible units where WWV is the WWV delay time, and S R is the range time measure- 
ment. There is an additional time correction applied in the AOPB differential correction program 
to the range time. When the range ambiguity is resolved the time is adjusted to correspond. 


A. 2 CORRECTION FOR THE AVERAGING EFFECT ON THE 7094 

After the smoothing process, a set of coefficients for the Chebyshev polynomial expressing 
the range rate measurement (§ RR ) as a function of raw data (station) time is available. 
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